/] /vlC r 

/ ~3 

/%~7 

77 


Final Report 

NASA Cooperative Agreement NCC2-450 
covering the period 
January, 1987 - January, 1989 
Principal Investigator - Jon A. Hoffmann 


The Influence of Free-Stream Turbulence 

on 

Turbulent Boundary Layers 
with 

Mild Adverse Pressure Gradients 


by 

J. A. Hoffmann, Professor 
S. M. Kassir, Graduate Student 
S. M. Larwood, Student Assistant 


November, 1988 


Aeronautical Engineering Department 
California Polytechnic State University 
San Luis Obispo, California 

<*ASA-CR- 184677) fill IHFlOIIiCI Cf N89-1S368 

fBEE-STB£AH T CBE OiBACI CM ICEECIi KT BQUMDABI 

iAXEES HUB RILE IDVEBSi PB£££0B£ GBADIEMIS 

final Report, Jan. 1967 - Jar. 1989 Unclas 

{California Polytechnic State Cciv.) 54 p G3/34 0187800 


\ 


Abstract 


The influence of near isotropic free -stream turbulence on the shape 
factors and skin friction coefficients of turbulent boundary layers is 
presented for the cases of zero and mild adverse pressure gradients. With 
free-stream turbulence, improved fluid mixing occurs in boundary layers with 
adverse pressure gradients relative to the zero pressure gradient condition, 
with the same free-stream turbulence intensity and length scale. Stronger 
boundary layers with lower shape factors occur as a result of a lower ratio of 
the integral scale of turbulence to the boundary layer thickness, and to 
vortex stretching of the turbulent eddies in the free-stream, both of which 
act to improve the transmission of momentum from the free-stream to the 
boundary layers . 
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Nomenclature 


A, B, D Constants for turbulence decay, dissipation length scale and 

autocorrelation coefficient equations 

Cf Boundary layer skin friction coefficient 

Acf, AH Change in Cf and H with free-stream turbulence compared to the 

low free-stream turbulence condition at the same a, Re^ and 
dP/dx. 


C 


P 


d 



H 

K 

LFST 

M 

n, m, k 
P 

IP- 

R* 

Re# 

t 

u' , v' , w' 
Uco 


Static pressure coefficient along plate, 

- P(x) -P(x-0 . 56m)/p (U w (x-0 . 56m) ) Z /2 

Diameter 

Time averaged dc voltage output of the hot wire anemometer system 

Time averaged rms voltage output of the hot wire anemometer 
system 

Boundary layer shape factor, - S*/9 
Constant, - Cf 0 (experimental)/Cf 0 (empirical) 

Dissipation length scale 

Low free-stream turbulence 

Mesh size of turbulence generators 

Exponents for turbulence decay, dissipation length scale and 
autocorrelation coefficient equations 

Static pressure 

T otal rms free-stream turbulence fluctuation, - 
/( u 72 + v 72 + v 72 )/! 

Free-stream Autocorrelation coefficient, - u'(t) u'(t -At)/ u'^ 

Reynolds number based on 6 , - Ha,9/v 

time 

Free-stream velocity fluctuations in the x, y and z directions 
respectively 

Free-stream time averaged velocity 
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a 

u* 


x,y,z 



a 


P 

*x 

6 , 6 *, $ 


0 

v 

r , 

P 

n 


Time averaged velocity 
Shear velocity 

Coordinate system relative to leading edge of flat plate 

Coordinate system relative to rods 

Distance from rods to leading edge of flat plate 

Angle of incidence of flat plate relative to the zero pressure 
gradient case. 

Pressure gradient parameter, - (6*/r 0 )/(d?/dx.) 

CO 

Integral length scale of turbulence, - J Rjj dx 

o 

Boundary layer thickness when a - 0.995 U,,,, displacement 
thickness and momentum thickness respectively. 

Divergence angle of diffuser wall 

Kinematic viscosity 

Wall shear stress 

Fluid density 

Wake parameter 


SUBSCRIPTS 

o Low free-stream turbulence condition 

MAX Maximum 
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Introduction 


Although many turbulent boundary layers encountered in engineering 
experience a turbulent flow above them, little information is available for 
designers which can be used to predict the effects that external turbulence 
will have on the boundary layers. Examples of flows with turbulent boundary 
layers under external turbulence or flows in which a boundary layer encounters 
the turbulent wake of a preceding object include turbomachinery blading, 
ejectors, airfoils downstream of a propeller or canard, as well as a wide 
variety of other internal and external flows where upstream free-stream 
turbulence exists. Also in nature, it is observed that fish swim in schools 
and birds fly in formation. 

Available literature has addressed the effects of isotropic homogeneous 
free-stream turbulence on the shape parameters of turbulent boundary layers 
with a zero pressure gradient. An investigation by Raghunathan (1) showed 
that the skin friction coefficient (cf) and the boundary layers thickness ( 6 ) 
increased, and the shape factor (H) decreased with increasing free-stream 
turbulence intensity cVPi/iu • Hancock (2) extended this result to 
incorporate the effect of the free-stream dissipation length scale (1^) , and 
found that Cf and S increased and H decreased when increasing a free-stream 
turbulence parameter defined as (Ju^/U^/dLx/S + 2). Castro (3) found that 
Hancock's results were Reynolds number dependent for Re^ < 2,000. 

The turbulent structure of a turbulent boundary layer, for the case of a 
zero pressure gradient, has been found to be affected by free-stream 
turbulence. An investigation by Huffman (4) for the case of isotropic free- 
stream turbulence found that the Reynolds stress and turbulent kinetic energy 
in the outer part of the boundary layer (y/5>0.4) increased as the free-stream 
turbulence intensity (J^/ U.) increased, and that these parameters were 
essentially independent of (vC^/iu for y/5<0.4. An investigation by Fang 
(5) for the case of near isotropic =1) and highly anisotropic free- 

stream turbulence J^/J^ > 2 .5) with a turbulent boundary layer in a near 
zero pressure gradient showed that the highly anisotropic free-stream 
turbulence condition results in larger isotropy (Vv^a/u 75 >1) in the outer 
part of the boundary layer (y /& > 0.6), and found that ~ 1 for the 

case of near isotropic free-stream turbulence when y/S > 0.6; for the case of 
a conventional turbulent boundary layer without free-stream turbulence. 
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< 1 for all y/S . The larger values of the Reynolds stresses, 
turbulent kinetic energy and isotropy ratios that exist in the outer part of 
boundary layers under free -stream turbulence, as reported by Huffman and Fang, 
are apparently factors which promote improved mixing and correspondingly 
decrease the shape factors of boundary layers relative to the low free- stream 
turbulence condition. 

The effects of free -stream turbulence on the shape parameters and the 
structure of turbulent boundary layers with adverse pressure gradients has 
been evaluated by Amal (17). Significant improvements in boundary layer 
strength were obtained when free -stream turbulence and adverse pressure 
gradients were coupled, resulting in a delay of separation which would 
otherwise occur without free-stream turbulence. 

The effects of near- isotropic and highly anisotropic inlet free-stream 
turbulence on the performance of a two-dimensional diffuser, an internal flow 
with an adverse pressure gradient, has been studied by Hoffmann (6) (7) (8). 
Improvements in the pressure recovery coefficient of the diffuser (Cp) of 
approximately 10% at a total included divergence angle (2©) of 12° (near peak 
Cp) and 24% at 2© - 20° (in the large transitory stall regime) were obtained. 
These improvements in Cp can be attained with a low turbulence generator total 
pressure loss; it was found that the ratio of the increase in Cp to the 
dimensionless total pressure loss of the turbulence generator could be larger 
than 13(8). These high gains were attained with high inlet total free- 

stream turbulence intensities (/q 72 /^ > 3.5%), high free-stream isotropy 
ratios > 2.1), and with the axes of the turbulence generators 

aligned so that the axes of the turbulent eddies are parallel to the diverging 
walls of the diffuser and perpendicular to the flow. It was also found that 
isotropic inlet free-stream turbulence can be significantly less effective at 
improving Cp as compared to highly anisotropic inlet free-stream turbulence 
(8). The increases in Cp were shown to be a result of both delayed separation 
and reduced distortion within the diffuser (7). 

In this study, the effects of free-stream turbulence on established 
turbulent boundary layers are evaluated. Free-stream turbulence can also 
alter the position of boundary layer transition; one recent paper on this 
topic is presented by Meier (14). 
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Generally, the most critical stage for a boundary layers development is 
encountered when it experiences a decelerating external free-stream flow and a 
corresponding adverse pressure gradient. For the case of boundary layers 
under low free-stream turbulence, this adverse pressure gradient can weaken 
the boundary layer (exhibited by an increase in the shape factor (H)) which 
may eventually lead to flow separation and a corresponding loss of pressure 
for internal flows, or a loss of lift and increase in drag for external flows. 
It is the purpose of this study to evaluate the effects of free-stream 
turbulence on turbulent boundary layers with mild adverse pressure gradients. 
The results will then be available for use by designers for both internal and 
external flows with free-stream turbulence. 
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Experimental System & Procedure 


To evaluate the effects of free-stream turbulence on turbulent boundary- 
layers, a flat plate was positioned in a wind tunnel and a set of rods was 
positioned upstream of the plate to generate free-stream turbulence. A 
detailed description of the system and instrumentation follows. 

All experiments were performed in the Cal Poly draw- thru wind tunnel, 
illustrated in Fig. 1. The test section has cross-sectional dimensions of 
0.88 mx 1.18 m. Screens and honeycomb upstream of the contraction, and a 
contraction ratio of 10:1, provide a relatively low turbulence intensity 
(/u^/U*, - 0.25X) at the test section. All experiments were performed 
holding the velocity upstream of the rods constant at 24.1 m/s. 

A set of monoplane parallel rods which span the entire cross-section of 
the test section of the wind tunnel with a solidity (d/M) of 0.32 and 
M - 6.67 cm. (see Fig. 2) was positioned at two locations upstream of the 
leading edge of the flat plate ( X Q - 0.67 m and X Q - 1.13 m) . 

The flat plate which spans the entire height of the test section of the 
wind tunnel as shown in Fig. 2 was constructed from 1.90 cm plywood, and was 
varnished and sanded to obtain a smooth surface. The plate pivots near the 
center (x - 0.63 m) to obtain various angles of incidence to the free-stream 
flow. The distance from the tunnel wall to the surface of the flat plate at 
the location of the pivot is 0.32 m. Experimentation was performed at three 
different angles of incidence. The plate was first positioned to obtain a 
zero static pressure gradient - the a -0* position; for a —0°, the plate 
and tunnel wall formed a diffuser with a 1.5* total included divergence angle. 
The plate was then adjusted to obtain adverse pressure gradients with 
a - 2.6° and a - 5.2°, inclined from the a - 0* position. The leading edge 
of the plate has a 2.5 cm radius; it was found that leading edge separation 
occurred at a - 2.6* and 5.2® when a 0.95 cm radius was used. Masking tape 
was applied to the leading edge of the plate, which tripped the boundary 
layer; leading edge separation also occurred at a — 2.6* and a — 5.2* when 
masking tape was not attached. 

The flat plate has six equally spaced static pressure taps located near 
the center line of the plate beginning at x - 0.30 m; at two of these x 
locations (x - 0.56m and x - 0.94 m) , the plate is instrumented with two 
Preston tubes of different diameter (d — 0.52mm and d — 0.82mm) and four fixed 
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total pressure probes with d - 0.82 mm located at various y positions used to 
determine velocities near the surface. The measured velocities near the 
surface fit the logarithmic law of the wall formula, equation 4, and were used 
to obtain u* and Cf for all test conditions. In order to determine velocities 
above the plate at any y position, a total pressure probe (d - 0.82 mm) was 
attached to a remote-controlled model B2509WIJ Velmex Inc. electric traversing 
mechanism. The voltage output of the mechanism, proportional to y, was 
measured using a HP model 3466A digital multimeter. The velocity pressure at 
each location was fed into a Celesco Model P15D± 0.1PSID pressure transducer, 
and then through a Model CD 10 Carrier Demodulator into a TSI model 1076 
multimeter, to obtain a ten second time average of the velocity pressure. A 
cubic spline curve fit of the velocity pressures and velocities was obtained, 
using a computer program, resulting in values of S, 6*, 0 and H for the 
boundary layers . 

Free -stream turbulence intensity and velocity traverses downstream of 
the rods were obtained using a TSI model 1050 anemometer and a TSI model 1052 
linearizer. The average voltage output of the linearizer (E) from a single 
wire normal to the mean flow, proportional to average velocity, was averaged 
for a 10 second period of time using a TSI model 1076 voltmeter, and the rms 
voltage output <y^> was averaged for a two second period of time using a HP 
model 3400A true rms meter. The probe was attached to a remote-controlled 
model B2509W1J Velmex Inc. electric traversing mechanism which produced a 
voltage output proportional to distance. These three voltages were fed into a 
HP model 1090A plotter to generate the velocity and turbulence vs. distance 
graphs in the non -homogenous flow region downstream of the rods. 

Free -stream turbulence measurements were obtained at two streamwise 
locations above the flat plate (x - 0.56 m and x - 0.94 m) using single wires 
normal to the flow and a wire inclined at 45* to the flow. The output of a 
single wire normal to the flow was also fed into a model 2032 B and K wave 
analyzer to generate plots of the autocorrelation coefficient vs. At. For 
these measurements, a TSI model 1057 signal conditioner was used after the 
linearizer, using a 50 kHz low pass filter to eliminate the noise of the 
anemometer and a 50 Hz high pass filter so that the values of more rapidly 
approach zero at large At, which thereby facilitated A x to be calculated more 
easily. The average decrease in the measured turbulence intensity 
due to the use of the high pass filter was 3.2%; no measurable difference 
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occurred due to the use of the low pass filter. Equations of the form 
Rx - De^ c , each valid for a given range of At, were obtained and these 
equations were then mathematically integrated to obtain the area under the 
curve and, correspondingly, the integral scale of turbulence. 


6 


Results and Discussion 


I) Free -Stream Flow Field Downstream of Rods 

The free -stream flow field downstream of the rods was first measured 
without the flat plate installed in the wind tunnel. For this condition, a 
small favorable pressure gradient exists due to boundary layer growth on the 
walls of the constant area test section. The average rms and dc voltage 
outputs of the linearizer of the hot wire anemometer system were obtained at 
several locations downstream of the rods, 0.38 < x'/M < 7.90, for 
0.27 ^ y'/M ^ 2.25 and are presented in Figs. 3 thru 8. The x'-y' axis is 
located at the center of the fourth rod from the tunnel wall (see Fig. 2). 

The average approach velocity for each traverse was 24.1 m/s. For x'/M < 7.9, 
the peak dc voltage output and the minimum rms voltage outputs of the hot 
wire, which occurs for flow between the rods, are directly proportional to u 
and /u' 2 respectively, whereas the voltage outputs at the other y' locations 
represent velocity and turbulence levels in the x'- y' plane. A small 
difference in the location of the peak velocity and the minimum turbulence 
intensity occurs for each traverse. This difference is attributable to the 10 
second time average used to obtain E and the two second time average used to 
obtain 

A plot representing the nonuniformity of the flow in the nonhomogeneous 
flow region downstream of the rods is presented in Fig. 9. Extrapolation of 
the plot shows that the variation of u as a function of y' is less than IX for 
x'/M > 8 or x' > 0.53m. In the region of essentially homogeneous flow 
downstream of the rods, turbulence decay was represented by equation 1; values 
of A and n are presented in table 1. 

u^/Ua, 2 - A x ' n Eqn. 1 

Now that the region with essentially homogeneous flow has been 
identified as x' > 0.53m, the flat plate was installed in the wind tunnel with 
its leading edge positioned in the homogeneous flow region downstream of the 
rods (X Q - 0.67m and 1.13m), and with a - 0® , 2.6* and 5.2® to obtain a zero 
pressure gradient and two conditions with mild adverse pressure gradients. 
Curves of the static pressure coefficient of the plate (this static pressure 
coefficient is based on the measured static pressure and free-stream velocity 
at x - 0.56 m) vs. x are presented in Fig. 10. The results show that the 
dimensionless static pressure gradient along the plate is essentially 
unaffected by the free-stream turbulence conditions and boundary layer growth 
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at each a, and is a weak function of x. A dimensionless pressure gradient 
parameter (/}) representing a ratio of the pressure forces to the shear forces 
in a section of the boundary layer (15) is presented in Table 2 and was not 
held constant in this study; dP/dx was calculated using a potential flow 
solution and the measured free -stream velocities at x - 0.56 m and x - 0.94 m. 
The turbulence measurements, presented in Table 2, show a near isotropic free- 
stream flow. For the case of a - 0* , the turbulence decay and the 
dissipation length scale equations (equations 1 and 2) were used to obtain the 
free-stream dissipation length scale (1^) as shown in equation 3. Values of 
A, B, n and m are presented in Table 1, and values of for each test 
condition at o — 0* are presented in Table 2. 

U a} d(u 72 )/dx' - -(u'’ 2 ) 3 / 2 /^ Eqn. 2 

Ljj - B x' m Eqn. 3 

The rate of turbulence decay - or the value of n in equation 1, has been 
found to be a function of the geometry, the solidity, and the distance 
downstream of the turbulence generator. For the zero pressure gradient 
condition, the values of A and n are listed in Table 1 and show that n 
increases as x' increases. 

Free-stream integral length scales (A x ) were obtained by finding the 
area under the Rjj vs. At curves; a sample curve from the wave analyzer 
presented in Fig. 11, shows a smooth curve representative of a broad-band 
turbulence. The results, presented in Table 2 and Fig. 12, show an increase 
in eddy size with downstream distance, and a decrease in eddy size with 
increasing a at any x' . The results indicate that vortex stretching occurs in 
the y direction due to diffusion in the y direction, which results in the 
lower measured values of A x at a - 2.6* and 5.2*, relative to a - 0* . 

II) Turbulent Boundary Layer Results 

The influence of free-stream turbulence on turbulent boundary layers is 
described below. The velocity profiles with different free-stream turbulence 
conditions and angles of plate incidence were the base boundary layer data 
obtained from which values of Cf, S and H were obtained. The results are 
presented in Table 2. As shown in Fig. 13, the velocity profiles near the 
wall fit the logarithmic zone law of the wall formula well (equation 4) , 
demonstrating that neither the free-stream turbulence levels nor the adverse 
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pressure gradients encountered in this study affect this relationship. The 
velocity profiles in the outer regions of the boundary layers have been 
related to a wake parameter, which has been shown to be dependent upon £ and a 
free-stream turbulence parameter (see the appendix). 

u /u* - 2.439 In u*y/V + 5.2 Eqn. 4 

A) Values of cj Q and H 0 

To evaluate the effects of free-stream turbulence on turbulent boundary 
layers, all values of cj and H for conditions with free-stream turbulence were 
compared to values of Cf Q and H 0 (values of Cf and H with low free-stream 
turbulence) holding Re ^ , a and dP/dx constant. Values of Cf and H, with and 
without free-stream turbulence, were compared at the same x locations at each 
a, where essentially the same values of dP/dx (or the same values of U,,, and 
dUjo/dx) exist. The data obtained resulted in different Re^ for boundary 
layers with and without free-stream turbulence at each x and a (as much as 
10X difference occurred) . Empirical formulas which fit the data of Cf Q and H Q 
at each value of x and a were used to obtain values of Cf Q and H 0 at the Re^ 
corresponding to each condition with free-stream turbulence. 

Values of H 0 at each a and x were obtained by numerical integration of 
formulas developed by von Doenhoff and Tetervin (9), equation 5 and 6. The 
division of equation 6 by equation 5, with U* and dU^dx held constant, 
presents a relationship for H 0 - f (0). The original work by von-Doenhoff and 
Tetervin ignored the effect of dP/dx on Cf. Equation 6, presented by Hirst 
(10) uses the Ludwieg and Tillman formula (11) (equation 7 with K-l) and 
thereby incorporates the effect of dP/dx on Cf. 

dff/dx - c fo / 2 - (H 0 +2)(tf/U <o ) dU.ydx Eqn. 5 

dH 0 /dx - e 4 - 68 < H o - 2.975) [ . (2/U,) (dU^dx) [5 . 89 log4.075Re tf ] 2 

-2.O35(H o -1.286)/0] Eqn. 6 

c fo - 0.246K / 10°- 678H oRe^ 0 - 268 Eqn. 7 

In the use of equations 5 and 6, the initial conditions used were the 
values of H 0 , d and Cf Q obtained at each x and a with low free-stream 
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turbulence. Values of dU^/dx were obtained from a potential flow solution 
using measured free-stream velocities at both x — 0.56 m and x - 0.94 m at 
each a. The maximum percentage change in H 0 due to this extrapolation, for 
the range of Re^ encountered in the tests with free-stream turbulence, was 
0.6X. 

Curves of H 0 - f(Re^) were also obtained from equation 5 and 6 , 
integrating with respect to x. The potential flow solution was used which 
presents U w and dU^dx - f(x). Curves which best fit the measured values of 
H Q and 9 at x - 0.56m and x - 0.94m at each a are presented in Fig. 14. All 
measured values of H Q fit these formulas within ± 1.5X. 

Values of Cf Q corresponding to the Re^ with free-stream turbulence for 
a — 0* at each x were obtained using a modified Clauser's formula (12), 
equation 8 , and values of cj 0 for o - 2.6* and 5.2* at each dP/dx were 
obtained using a modified Ludwieg and Tillman formula (11), equation 7. In 
these equations, values of H 0 were obtained from the von Doenhoff and Tetervin 
formulas . 

c fo - K/(3 . 96 log H o Re 0 + 3. 04) 2 Eqn. 8 

The value of K is an indicator of the magnitude of the departure of the 
results of this study from results with K - 1 as obtained by Clauser, and 
Ludwieg and Tillman. As shown in Table 2a, the maximum departure obtained was 
5X. The maximum precentage change in Cf Q due to this extrapolation, for the 
range of Re^ encountered in the tests with free-stream turbulence, was 1.3X. 

Curves of Cf D — f (Re^) which best fit the data at both x - 0.56m and 

x — 0.94m were obtained using equation 8 for a — 0* and equation 7 for 

a - 2.6* and 5.2*. The results are presented in Fig. 15. All measured values 

of Cf Q fit these formulas within ± 1.6X. 

The magnitudes of H Q and Cf Q at each a, as shown in Figs. 14 and 15, are 
affected by different initial conditions for the boundary layers which occur 
at each a. It is observed that the value of H 0 at x - 0.56m and a - 5.2* is 
lower compared to the value at x - 0.56m and at - 2.6*, and the trend of Cf 0 at 
a - 5.2* appears high relative to that at a - 2.6*. Larger values of U,,, over 
the plate occur as a increases due to blockage caused by the plate in the test 
section; it is concluded that this flow acceleration results in lower values 
of H 0 and larger values of Cf Q near the leading edge of the plate with 
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increasing a, as compared to the values that would occur without plate 
blockage. In this study, the boundary layers with free-stream turbulence are 
compared to those without free-stream turbulence at the same values of a (as 
well as x and Re^) and therefore with essentially the same free-stream 
velocity distribution. 

B) Zero Pressure Gradient Case with Free-Stream Turbulence. 

For the conditions with free-stream turbulence and a zero pressure 
gradient, the results of this study are compared to those of Hancock (2). A 
comparison of values of Acf/cf Q at various values of (L^/S + 2) is 

presented in Fig. 16 and the relationship between Acf/cf Q and AH/Ho is 
presented in Fig. 17. The results of this study agree reasonably well with 
those of Hancock (2), considering differences in the experimental system used 
by both investigators (e.g. Hancock used a flat plate with a sharp leading 
edge, while a rounded leading edge with some flow acceleration after the rods 
due to plate blockage was used in this study) and the repeatability of the 
results of both investigations. 

C) Adverse Pressure Gradient Case with Free-Stream Turbulence. 

Curves of H vs. Re^ for a - 0° , 2.6° and 5.2° are presented in Fig. 14. 
For the case with free-stream turbulence, relatively low values of H occur for 
a - 2.6° and a - 5.2° compared to a - 0° for constant J q^/u«; values of H at 
a - 2.6° and a - 5.2° appear to be relatively independent of Re^ for constant 
Traditionally for the case without free-stream turbulence, an 
increase in H is expected with increasing Re^, and the rate of increase of H 
should increase as a increases, as shown on Fig. 14 for the cases without 
free-stream turbulence. Curves of Cf vs. Re^ shown in Fig. 15 also show large 
magnitudes of Cf for conditions with free-stream turbulence at all a, relative 
to the values of Cf Q . An explanation for these remarkable improvements in 
boundary layer strength is addressed below. 

In Fig. 18, values of AH/H Q are plotted as a function of a free-stream 
turbulence parameter, , based on the integral length scale. This 

parameter, which places a larger weight on the free-stream length scale, is 
somewhat different than the free-stream turbulence parameter used by Hancock, 
Tu^/U^ /(L^/S + 2), which is based on the dissipation length scale. In Fig. 
18, it can be seen that a best fit curve for a - 0° exists, and a separate 
best fit curve for both a - 2.6° and a - 5.2° exists. As shown in Fig. 19, 
the data for Acf/cf 0 is again generally lower for a - 0° , relative to the data 
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for both a - 2.6° and a - 5.2°. With 7q^ - 2.1, values of AH/H Q and 
Acf/Cf 0 at a - 2.6° are both about 25 % larger compared to those at a - 0° . 

As a result of a higher repeatability of Acf/cf Q relative to AH/H Q , more 
distinction of the results of this study can be obtained by evaluating the 
data of Acf/cf 0 . An inspection of Fig. 19 shows that at a - 0* , values of 
Acf/cf 0 appear to be independent of Re^ for 2500<Re^<3900. Castro (3) found 
that for a - 0* , values of Acf/cf 0 were independent of Re^ for Re^>2000. The 
single curve presented for a - 2.6° and a - 5.2° appears to be applicable for 
Re0“5OOO when a - 2.6° and for Re^“9000 when a - 5.2°; the results at lower 
Re^ with a — 5.2" fall below this curve. More data with a larger range of Re^ 
needs to be obtained to better define the limitations of this preliminary 
study. 

The results presented suggest that an improvement in fluid mixing occurs 
for both a - 2.6* and a - 5.2*, relative to a - 0° . Two major differences 
occur with boundary layers under free-stream turbulence with adverse pressure 
gradients relative to the zero pressure gradient condition. First, the 
boundary layer thickness is larger for the adverse pressure gradient condition 
relative to the zero pressure gradient condition. For the same free-stream 
length scale, values of X x /S are smaller resulting in more effective mixing. 
Hancock (2) has also shown that lower values of I^/S result in larger values 
of AH/H q and ACf/Cf Q for the condition of a - 0° with the same Tq^/U^. Since 
the free-stream turbulence parameter reflects changes in both Jq* 2 /U 00 and 
A x /$, one would expect that another fluid mechanism exists which causes the 
larger values of AH/H Q and Acf/cf 0 at constant Jq^/V^Xx/S with adverse 
pressure gradients relative to the zero pressure gradient condition. Second, 
for the case of boundary layers with an adverse pressure gradient, vortex 
stretching of the turbulent eddies in the free-stream occurs in the y 
direction as a result of diffusion of the flow in the y direction; this is 
reflected by a decrease in A x with increasing a as shown in Fig. 12. It is 
concluded that the effectiveness of fluid mixing is augmented for conditions 
with adverse pressure gradients and with free-stream turbulence, due to vortex 
stretching and an associated increase in vorticity of the eddies. In another 
study, Viets (13) investigated the effect of vortex stretching within the 
boundary layer for flow in a diffuser with low free-stream inlet turbulence. 

He found that diffusion normal to the plane of the velocity profile leads to 
improved mixing in contrast to diffusion in the plane of the velocity profile. 
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For this study, improved mixing occurs with free-stream turbulence and 
diffusion in the plane of the velocity profile. The work of Arnal (17) has 
also revealed that increases in boundary layer strength can occur when free- 
stream turbulence and adverse pressure gradients are coupled; the work 
presented in this paper identifies, for the first time, the fluid mechanisms 
responsible for this improvement in boundary layer strength. In Arnal' s 
experiment with a large adverse pressure gradient and a boundary layer near 
separation for the low free-stream turbulence condition, relatively large 
fractional changes in H and cj were obtained AH/H Q “0.25 and ACf/Cf 0 m 5) 
when was approximately 5%. A direct comparison of the results of 

Arnal with the results presented in this paper is made difficult due to the 
lack of A x data in Arnal' s work. 

A plot of Acf/Cf 0 vs. AH/H 0 , presented in Fig. 17, shows that for the 
same incremental change in Acf/cf 0 lower values of AH/H 0 generally occur for a 
-2.6* and a - 5.2*, relative to Hancock's results for a - 0* . The two data 
points on the left of Hancock's curve occur for low Re^ for the respective 
values of a; Castro (3) also obtained data on the left of Hancock's curve for 
low Re0 with a - 0* . Castro (3) also obtained data on the left of Hancock's 
curve for low Re^ with a - 0* . 

Ill) Accuracy and Repeatability 

In this study, values of H and Cf are related to a free-stream 
turbulence parameter, 7q' ^/U^/A jj/ 5 . Values of H and Cf were both found to be 
repeatable within ± 1.5%, with 20:1 odds. As a result of the relatively small 
values of Acf and AH, the repeatability of Acf/Cf 0 is ± 0.02, or ± 10% with 
ACf/Cf 0 - 0.2. The repeatability of AH/H 0 is ± 0.02, or ± 20% if AH/H Q is 
0.1. Actual magnitudes of Cf for the condition of low free-stream turbulence 
were found to differ from those of other investigation by as much as 5%, 
probably due to differences in experimental systems. Values of the free- 
stream turbulence parameter were found to be repeatable within ± 10%. 
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Conclusions 


Traditionally, adverse pressure gradients have been considered a menace 
to the boundary layer which can weaken it and eventually lead to flow 
separation. Indeed, this is true for the case of a boundary layer with low 
free-stream turbulence. The results of this study have shown that, for the 
case of turbulent boundary layers under external free-stream turbulence, 
improved mixing and a more effective transmission of momentum from the free- 
stream to the boundary layer occurs when mild adverse pressure gradients 
exist, relative to the zero pressure gradient case with the same values of the 
free-stream total turbulence intensity and integral length scale. The 
boundary layer shape factor has been shown to remain constant, within the 
limits of the uncertainity of the data, with increasing Re^. The fluid 
mechanisms responsible for this improved mixing and corresponding stronger 
boundary layers with lower shape factors are: 

1) A reduction of the ratio of the free-stream integral scale of 
turbulence to the boundary layer thickness (A x /$) generally occurs with 
boundary layers in adverse pressure gradient flow fields. This 
reduction in A x /$ is primarily due to the increase in S which naturally 
occurs when adverse pressure gradients exist, relative to smaller values 
of 6 which occur for the zero pressure gradient case. 

2) Vortex stretching of the turbulent eddies in the free-stream 
occurs when adverse pressure gradients exist, resulting in a higher 
vorticity and a corresponding augmentation of fluid mixing. 

In this manuscript, values of the fractional changes in H and cj that 
occur with free-stream turbulence relative to the values without free-stream 
turbulence, are presented as a function of a free-stream turbulence parameter, 
Jq 9 2 /U 00 /A x /S , for conditions with zero and mild adverse pressure gradients. 

The changes in cj and H obtained appear to be dependent upon Re^ for low Re^. 
For this preliminary study, values of A x /$ were on the order of magnitude of 
1. Intuitively, one would expect that the free-stream would not exert as much 
effect on the boundary layer for very small values if A x /$ ; for very large 
values of A x /5 , an unsteady boundary layer would be expected. Future studies 
should investigate the effect of free-stream turbulence on boundary layers 
with favorable pressure gradients, as well as boundary layers with high 
adverse pressure gradients, using a larger range of free-stream turbulence 
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conditions as well as a larger range of Re^ . Now that the fluid mechanisms 
resulting in improved boundary layer strength have been identified, 
mathematical modeling can begin to predict boundary layer parameters for a 
wide variety of flow conditions with free-stream turbulence. 
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Table 1 


Constants and Coefficients for Turbulence Decay 

and 

Integral Length Scale Equations 


Condition A n B m 

Test Section w/o Plate 


0.50m < 

x' < 1.27m 

16.34 

-1.79 

2.26 

0.105 

a - 0°, 

X Q - 0.67m 





1.23m < 

x'< 1.61m 

3.87 

-1.542 

1.276 

0.229 

a - O', 

X Q - 1 . 13m 





1 . 69m < 

x' < 2.07m 

0.125 

-0.897 

0.394 

0.552 


NOTE: For equations 1 and 3, x' and A x have the units of cm. 
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Table 2a 


Summary of Results 


Condition 

X 

(m) 

u» 

(m/s) 

S 

(cm) 

9 

(cm) 

Ho 

C f 

•‘■O 

xlO 3 

K 

Re# 

P 


0.56 

27.2 

1.47 

0.156 

1.410 

3.33 

0.998 

2850 

0 

a - 0° 










LFST 

0.94 

27.0 

2.01 

0.229 

1.420 

2.99 

0.965 

4152 

0 


0.56 

29.9 

1.68 

0.201 

1.455 

2.87 

1.047 

4038 

0.33 

cr - 2.6° 










LFST 

0.94 

28.4 

2.77 

0.338 

1.443 

2.44 

0.990 

6452 

0.42 


0.56 

32.6 

2.36 

0.297 

1.431 

2.63 

1.051 

6521 

0.79 

a - 5.2° 










LFST 

0.94 

30.3 

3.45 

0.462 

1.480 

2.12 

1.009 

9418 

1.41 


20 


Table 2b 


Summary of Results 











Jn™/ u. 

(Lx/7+7) 

Condition 

X 

(ffl) 

<%) 



u» 

(%) 

Ueo 

(m/s) 

^x 

(cm) 

A x 

(cm) 

a - 0° 

0.56 

4.83 

1.10 

1.27 

5.45 

26.3 

3.84 

5.48 

0.0133 

X D - 0.67m 

0.94 

3.92 

1.17 

1.38 

4.68 

26.1 

4.09 

6.09 

0.0117 

a - 2.6° 

0.56 

4.67 

1.09 

1.05 

4.88 

29.2 


5.37 


X Q - 0.67m 

0.94 

4.23 

0.86 

0.90 

3.90 

27.6 


5.92 


o - 5 . 2* 

0.56 

4.63 

1.23 

1.06 

5.10 

32.4 


5.15 


X Q - 0.67m 

0.94 

4.16 

1.03 

1.02 

4.29 

29.5 


5.84 


a - 0° 

0.56 

3.55 

1.10 

1.27 

4.00 

25.7 

6.69 

6.20 

0.0072 

X c - 1.13m 

0.94 

3.24 

0.94 

1.16 

3.37 

25.6 

7.49 

6.47 

0.0071 

a - 2.6* 

0.56 

3.38 

1.32 

1.14 

3.92 

29.7 


5.96 


X Q - 1.13m 

0.94 

3.11 

1.19 

1.13 

3.45 

27.7 


6.22 


a - 5.2° 

0.56 

3.19 

1.35 

1.18 

3.78 

32.5 


5.98 


X Q - 1.13m 

0.94 

3.13 

1.28 

1.10 

3.55 

29.8 


6.05 
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Table 2c 


Summary of Results 


Condition 

X 

(m) 

S 

(cm) 

9 

(cm) 

H 

c f 

xlO 3 

V 

A'V 

-TJT- 

(%) 

_ AH 
H o 

Acf 


0.56 

2.36 

0.146 

1.326 

3.93 

2582 

2.35 

0.067 

0.164 

a - 0° 










X„ - 0.67m 

0.94 

3.00 

0.219 

1.322 

3.50 

3972 

2.31 

0.072 

0.164 


0.56 

3.15 

0.208 

1.318 

3.52 

4082 

2.86 

0.094 

0.230 

a - 2.6° 










X, - 0.67m 

0.94 

4.90 

0.333 

1.302 

3.13 

6185 

3.23 

0.096 

0.261 


0.56 

3.75 

0.298 

1.302 

3.19 

6499 

3.71 

0.090 

0.211 

a - 5.2* 










X„ - 0.67m 

0.94 

5.50 

0.458 

1.305 

2.80 

9108 

4.04 

0.114 

0.295 


0.56 

2.26 

0.155 

1.347 

3.73 

2718 

1.46 

0.047 

0.112 

a - 0° 










X„ - 1.13m 

0.94 

2.95 

0.224 

1.343 

3.39 

3978 

1.53 

0.057 

0.127 


0.56 

3.23 

0.211 

1.310 

3.38 

4207 

2.12 

0.100 

0.191 

a - 2.6° 










X, - 1.13m 

0.94 

4.58 

0.348 

1.321 

2.97 

6478 

2.54 

0.085 

0.219 


0.56 

3.82 

0.292 

1.311 

2.99 

6380 

2.41 

0.083 

0.125 

a - 5.2* 










X„ - 1.13m 

0.94 

5.55 

0.467 

1.324 

2.64 

9352 

3.26 

0.105 

0.240 
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Table 2d 


Summary of Results 


Condition 

X 

(m) 

c fo 

xlO 3 

H 0 

P 

a - 0° 

0.56 

3.377 

1.418 

0 

X D - 0.67m 

0.94 

3.008 

1.424 

0 

a - 2.6* 

0.56 

2.861 

1.455 

0.26 

X Q - 0.67m 

0.94 

2.482 

1.440 

0.43 

a - 5.2* 

0.56 

2.634 

1.431 

0.60 

X Q - 0.67m 

0.94 

2.162 

1.474 

0.94 

a - 0* 

0.56 

3.353 

1.414 

0 

X Q - 1.13m 

0.94 

3.007 

1.424 

0 

a - 2.6* 

0.56 

2.837 

1.455 

0.27 

X Q - 1.13m 

0.94 

2.436 

1.443 

0.48 

a - 5.2° 

0.56 

2.657 

1.429 

0.65 

X Q - 1.13m 

0.94 

2.129 

1.479 

1.03 
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APPENDIX 


Velocity Profile Correlation 


The velocity profiles of the boundary layers were found to fit a formula 
suggested by Hancock (2) and Dean (21), equation 9. 

a/u* - 2.439 In u* y/v + 5.2 + 2.439[(1 + 6H)(y/5) 2 - (1 + 4H)(y/5) 3 ] 

Eqn. 9 

The first two terms of eqn. 9 represent the law of the wall, equation 4, 
whereas the last term represents the deviation from the law of the wall (the 
wake function) and becomes significant in the wake region of the boundary 
layers. In this appendix, values of II , the wake parameter, are correlated 
with fi and free-stream turbulence parameters. The velocity profiles obtained 
were found to fit eqn. 9 reasonable well for both positive and negative values 
of II, whereas they did not fit Coles' wake function (2 sin 2 (tt/ 2 ) (y/5)) well 
for negative values of II. 

a) For the case of fi - 0 and low free-stream turbulence, the values 
of II fit the results presented by Cebeci (15) reasonably well, as 
shown in Fig. 20. 

b) For the case of fi - 0 and with free-stream turbulence, the results 
of Hancock (2) and the present results, presented in Fig. 21, show 
that H is a st rong function of the free-stream turbulence 
parameter ,fu' ^ /I!*, /((I^ / 5) +2). 

c) For the case of fi > 0, values of n are shown to be a strong 
function of fi (see Fig. 22). The results of this study with low 
free-stream turbulence fit the results presented by White (19) II - 
0.8 (0.5 + £)0* 23 and Das (20) for equilibrium boundary layers 
reasonably well. 

d) For the case of fi > 0 and with free-stream turbulence, the results 
obtained by Evans (16) along with the present results form a 
distinct separate curve which parallels the curves obtained by 
White (19) and Das (20). The results of Arnal (17), both with and 
without free-stream turbulence, do not fit any of the curves 
presented. These results then suggest that II is a function of 
both fi and a free-stream turbulence parameter. A reasonably good 
correlation of these parameters is presented in Fig. 23. More 
data needs to be obtained to better define the curves in this 
graph, and to obtain this correlation for larger values of fi. 
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Fig. 2 

Flat Plate and Turbulence Generator 
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FIG. 10 

Plate static pressure coefficient along plate 
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0-5.2°, X Q - 1.13m, x - 0.56m 
Fig. 11 

Autocorrelation Coefficient vs. At 
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FIG. 13 

Velocity profiles in logarithmic law of the wall region 






FIG. 14 

Boundary layer shape factor vs. Reynolds number 
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FIG. 15 

Boundary layer skin friction coefficient vs. Reynolds number 






FIG. 16 

Fractional change of skin friction vs. Hancock’s free-stream turbulence parameter 
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FIG. 17 

Fractional change in skin friction vs. fractional change in shape factor 
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FIG. 19 


Fractional change in skin friction vs. free-stream turbulence parameter 
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Fig. 20 Variation of wake parameter 
with momentum thickness Reynolds number 
for 6=0° and LFST. 





Hancock (2) (17OO<Re0<6OOO) 
Present Results 
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Fig. 21 Wake parameter as a function of a free-stream turbulence 
parameter for p = 0 






Fig. 23 Wake paramoter as a function of turbulence intensity for constant p 




